Substrate-induced bianisotropy in metamaterials 
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We demonstrate that the presence of a supporting substrate can break the symmetry of a metamaterial 
structure, changing the symmetry of its effective parameters, and giving rise to bianisotropy. This indicates 
that magneto-electric coupling will occur in all metamaterials fabricated on a substrate, including those with 
symmetric designs. 



Optical metamaterials usually consist of a patterned 
metal-dielectric composite structure mounted on a sub- 
strate. They are designed to exhibit exotic properties, 
such as a negative refractive indejP'. The design pro- 
cess usually aims to achieve the maximum possible sym- 
metry, to better approximate an ideal isotropic mate- 
rial. In fabricating optical metamaterials, it is impor- 
tant to have straight side walls, since tapering breaks 
symmetry and causes bianisotropy. A bianisotropic ma- 
terial is one which develops an electric polarization in 
response to a magnetic field, and vice-versa, with mate- 
rial parameters of the form Di = eoCijEj — i^ijHj/c and 
Bi = noHijHj +i£^fjEj/c. This complicates the effective 
medium description, and can inhibit the negative refrac- 
tive index"^. The effect of fabrication imperfections has 
been studied in order to characterise bianisotropy with 
the aim to reduce its effecl!^. However, in this Letter we 
show that the presence of a substrate intrinsically induces 
bianisotropy in a metamaterial, which should be taken 
into account for accurate characterization and control of 
the metamaterial properties. 

A substrate has previously been shown to significantly 
influence the plasmonic resonances associated with a neg- 
ative refractive inde>!P. However, to date little attention 
has been paid to the resultant changes in symmetry, an 
important exception being "planar chiral" structured. 
In this case an essentially two-dimensional structure can 
exhibit optical activity, either when the metamaterial ar- 
ray has low in-plane symmetry, or when the sample is not 
normal to the incident wave- vector. We will show here a 
different form of substrate-induced symmetry breaking. 
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FIG. 1. (a) System under consideration, and (b) representa- 
tion as hybridization of non-identical resonators 
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for three-dimensional symmetric structures illuminated 
at normal incidence. 

The negative index in optical metamaterials arises 
when there is both a negative electric and magnetic re- 
sponse. The negative electric response occurs naturally 
in metals, while the negative magnetic response requires 
a pair of metal structures having an anti-symmetric plas- 
monic mode. Thus we consider here a general model of 
a pair of plasmonic resonators, which could be the two 
metal layers of a fishnet, a cut-wire pair or a plasmonic 
dimer. In Fig. [ija) we show the structure under con- 
sideration, noting that the details of this patterning are 
unimportant for this argument, except that the pattern 
is uniform through the layers. 

Conceptually, bringing the two identical patterned 
metallic layers together will cause their resonances to hy- 
bridize into symmetric and anti-symmetric modes. These 
will then respond to the imposed electric and magnetic 
field respectively. In practice there is substantial cou- 
pling of evanescent waves into the dielectric substrate, 
and this coupling will be different for the two plasmonic 
layers. Therefore we can consider that the system is 
formed by the hybridization of two non-identical plas- 
monic resonators, as shown in Fig. [ijb). The nominally 
symmetric mode will have some anti-symmetric compo- 
nent, and vice- versa. Each mode will therefore couple 
to both electric and magnetic fields, resulting in bian- 
isotropy. The argument is essentially identical to that 
for intrinsically non-symmetric structures, once it is ap- 
preciated that the substrate breaks the symmetry of the 
system, even if it is outside the boundaries imposed dur- 
ing the retrieval procedure. We will show a specific ex- 
ample of the experimentally important fishnet structure. 
We consider the case of a semi-infinite substrate, however 
the results still apply for a finite substrate. 

The material parameters of a homogeneous structure 
can be retrieved from transmission and reflection datcP. 
By utilizing the two values of reflection from different 
sides of the structure, this approach has been extended to 
asymmetric structures exhibiting bianisotropy^"—. Af- 
ter calculating the generalized scattering parameters, we 
consider the two values of the reflection coefficient when 
illuminating from free space (5*11) or through the sub- 
strate (S'22). Even for a symmetric structure, these 
will not be equal, since they are defined with respect 
to dilferent values of the wave impedance. However, 
we can transform the scattering matrix S with refer- 
ence impedances (?7oi'7s«6) to S" with both reference 
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(c) Extracted parameters 
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FIG. 2. Simulation results of the fishnet structure (a) Trans- 
mission and reflection magnitude (b) real part of extracted 
impedance values (c) real part of extracted permittivity, per- 
meability, and magneto-electric coupling. The dashed lines 
indicate the (approximately) anti-symmetric plasmonic reso- 
nances. 



impedances being 779 , and we find that for a simple dielec- 
tric layer, the reflection coefficients become equal, and 
this should also be true for a symmetric metamaterial. 

We then apply this approach to a fishnet struc- 
ture, with transverse period of 400nm, consisting of 
200x350nm holes in metal-dielectric-metal layers with 
metal thickness 30nm and dielectric thickness 30nm, on 
a substrate of index 1.5 which is impedance matched to 
free space (chosen for consistency with subsequent re- 
sults, as justified below). The structure was modeled 
in CST Microwave StudicJi^, with the metal being mod- 
eled using a 4th order fit to the parameters of silver, 
and the dielectric layer having e = 2.75. The transmis- 
sion of this structure is shown in Fig. [2][a) to exhibit two 
dips, corresponding to the excitation of approximately 
anti-symmetric plasmonic modes, and a maximum corre- 
sponding to the cut-off resonance of the hole mod^^^. ft 
can be clearly seen that the two reflection coefficients are 



different in magnitude, with the most significant differ- 
ences corresponding to the frequencies where plasmonic 
modes are excited. We note that in modeling of lossless 
structures, bianisotropy may not be immediately appar- 
ent, since 15*111 = |'5'22| is enforced by energy conservation 
and reciprocity. However, the phase of these two param- 
eters may be altered by the substrate, hence bianisotropy 
can still occur. 

We apply the method given in Ref.|TT]to extract the the 
two values of impedance, which are plotted in Fig. [2]^b). 
Failing to account for the bianisotropy in the extraction 
results in some error in the effective index, and yields 
an impedance value close to that of z\ (not shown). It 
is clear that the variation in impedance is strongest near 
the plasmonic features, supporting the argument that the 
non-symmetric hybridization of plasmonic modes is the 
cause of bianisotropy. In Fig. |2][c) we extract the equiv- 
alent parameters of the structure using the method from 
Ref. [HI where for clarity only the real parts are shown. It 
can be seen that the relevant component of the magneto- 
electric coupling, ^a;j, is significant around the magnetic 
resonance. The observed values indicate that the mag- 
netic polarization excited by the electric field is about a 
third as strong as that excited by the magnetic field. In- 
terestingly, at ^260TIIz, which corresponds to the cut-off 
resonance of the holes where e there is a strong vari- 
ation between impedance values, but no corresponding 
feature in the extracted magneto-electric coupling. 

In Fig. [3] we plot the maximum magnitude of i^xy 
over the considered frequency range, as a function of 
the substrate parameters. Whilst in optical experiments 
the substrate would have a purely dielectric response, 
we consider also a substrate with a magnetic response 
to better understand the substrate influence. Increas- 
ing tsuh or /ijub results in higher bianisotropy, although 
tsuh has a much greater influence. The dashed line in- 
dicates the impedance match between the substrate and 
free space, and the absence of any features in this re- 
gion indicates that bianisotropy cannot be attributed to 
an impedance matching effect. On the other hand, for 
flxed impedance or permeability bianisotropy increases 
monotonically with the index of the substrate. This is 
consistent with our argument that the change in plas- 
mon dispersion is the cause of this effect. 

While the extracted parameters reproduce the original 
transmission and reflection data up to numerical accu- 
racy, they are non-local in nature^. Therefore we shall 
probe the internal field of a structure directly, in order to 
demonstrate that the magneto-electric polarization is a 
real physical effect which does not rely on the homogeniz- 
ability of the structure. With reference to Fig. [ija), we 
excite the structure simultaneously from both directions 
with normally incident waves. This creates a standing 
wave pattern, and by adjusting the relative phase of the 
incident waves it is possible to position the center of the 
structure at a zero of the electric or magnetic field. Given 
that the thickness of the structure under consideration is 
much less than the wavelength, we have a good approx- 
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FIG. 3. Maximum real part of ^^,1, as a function of sub- 
strate permittivity and permeability. The dashed line indi- 
cates impedance matching to free space. 
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FIG. 4. Magnitude of the magnetic dipole moment induced in 
a fishnet under magnetic (rrim) and electric (me) excitation, 
with dashed lines showing the plasmon resonance frequencies. 



imation of pure excitation by the magnetic or electric 
field respectively. In order to simplify the excitation of 
the standing wave and the interpretation of the results, 
we consider the case where the substrate impedance is 
matched to free space, and the index is 1.5. We calcu- 
late the induced magnetic and electric dipole moments by 
evaluating the integrals m — juj/2 J x x {e — l)E(Px and 
p = J{e — l)E(Px over the metal and dielectric regions. 

The X component of the induced magnetic dipole mo- 
ment is shown in Fig. |4] for excitation by an electric field 
TOg, and is compared to rrim: the magnetic polarization 
induced by the magnetic field. Quantities are normal- 
ized to the volume of the unit cell, and to the incident 
magnetic field, thus they are dimensionally equivalent to 
polarisabilities Xme and Xmm- It is clearly observable 
that there are maxima of magneto-electric coupling cor- 
responding to the frequencies plasmon excitation observ- 
able in Fig.jga). Si gnificantly, at some frequencies, the 
magneto-electric polarization is stronger than the purely 
magnetic polarization. In addition, there is also some 
background cross-polarization which increases with fre- 



quency, which is not observable in the effective parame- 
ters. 

We expect that the effect predicted and analyzed here 
should occur in multi-layered structures, however in this 
case the analysis becomes more complicated. In particu- 
lar, since this is an interface effect, the multi-valued na- 
ture of the impedance should not strongly depend on the 
length of the structure. However, inversion of the trans- 
mission and refiection assumes homogeneous parameters, 
so the effect will be incorrectly attributed to bianisotropy 
throughout the bulk. We note that approaches which ac- 
count for surface layers of a metamaterial exislfSi^ but 
are not widely utilized in the literature. We further pro- 
pose that at the interface between a bulk metamaterial 
and free-space, the broken symmetry could also result 
in local bianisotropy. This would not be observable by 
any technique based on S-parameter inversion of a finite 
thickness slab, since the effect on both reflection coeffi- 
cients would be identical. 

In conclusion, we have demonstrated that the substrate 
breaks the symmetry of structurally-symmetric metama- 
terials, resulting in significant bianisotropy. We have 
shown that this effect is strongest at the resonance of 
plasmonic modes, which can be understood as the hy- 
bridization of non identical resonators to form the meta- 
material. Our results suggest that most metamaterials 
reported in the literature will exhibit bianisotropy, and 
should be taken into account for accurate analysis of the 
metamaterial properties. 
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